Unconventional Magnetism in a Nitrogen-Containing Analogue of Cupric Oxide 
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We have investigated the magnetic properties of CuNCN, the first nitrogen-based analogue of 
cupric oxide, CuO. Our muon spin relaxation, nuclear magnetic resonance and electron spin reso- 
nance studies reveal that classical magnetic ordering is absent down to lowest temperatures. How- 
ever, large enhancement of spin correlations and unexpected inhomogeneous magnetism have been 
observed below 80 K. We attribute this to a peculiar fragility of the electronic state against weak 
perturbations due to geometrical frustration, which selects between competing spin-liquid and more 
conventional frozen states. 

PACS numbers: 75.10.Kt, 76.75.-K, 76.60.-k, 76.30.-v 



Numerous copper oxides fascinate with a plethora 
of exotic electronic ground states, including high- 
temperature superconductivity and quantum magnetism 
[3, El. Many of them arc low-dimensional spin- 1/2 geo- 
metrically frustrated antiferromagnets (AFM) Q. There, 
quantum fluctuations compete with conventional long- 
range ordering (LRO) and can prevail down to zero tem- 
perature, leading to cooperative quantum ground states 
with spin-liquid (SL) character. In copper oxides, a SL 
ground state has recently been experimentally reported 
on two-dimensional (2D) kagome and triangular lattices 
d, [B| . In the latter case, numerous competing phases - 
various SL states, (in)commcnsurate LRO - have been 
theoretically predicted for isosceles triangles, depending 
on the ratio of the two exchange constants [6H8| . 

Replacing oxygen with other anions is an important 
step in tailoring their physical properties. Transition- 
metal carbodiimides, M(NCN), are good examples for 
an "organic" route. They are analogues of oxides because 
the complex anion NCN 2 " is isolobal with O 2 - @. Their 
existence was elusive until very recently, when metathesis 
reactions lead to the entire family, M = Mn, Fe, Co, Ni, 
and Cu [13, El- The crystal structure of CuNCN (Fig. [J) 
is orthorhombic (Cmcm) and consists of corrugated lay- 
ers with a 4+2 nitrogen environment of Cu 2+ [l2|, [H| . 
Unlike other M(NCN) materials, all being AFMs with 
a Neel temperature around 200-300 K, CuNCN exhibits 
a nearly temperature (T) independent and surprisingly 
small (xb ~ 9 • 10 -5 emu/mol) bulk magnetic susceptibil- 
ity, which enhances at low T, presumably due to a tiny 
impurity amount 



ner predicted a quasi-one-dimensional magnetism, with 
the leading AFM interaction J' ps 2500 K bridged via 
NCN [see Fig. QTa)] and LRO around 100 K, due to 
a ferromagnetic interchain interaction J\ [l3|. Neither 
bulk magnetic nor neutron diffraction (including polar- 
ized 13J) measurements, however, could confirm LRO. 
Liu et al. [l2| had proposed an alternative model, treat- 
ing CuNCN as the 2D spatially anisotropic triangular 
lattice in the ab plane [Fig.QJb)]. Their GGA+C7 calcu- 
lations again predicted strong AFM exchange couplings 
(Ji 800-1000 K) mediated by N atoms. Their size 
and the frustrated geometry then lead Tchougreeff and 
Dronskowski to propose a 2D resonating-valence-bond 
(RVB) picture with a nonmagnetic SL ground state [Hj]. 
These ambiguities about the ground state call for detailed 
microscopic investigations, allowing detection of intrin- 
sic spin-only susceptibility. In this Letter we present 
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The origin of the surprising reduction in Xb has been 
debated. Based on LDA calculations, Tsirlin and Ros- 



FIG. 1: (color online), (a) Unit cell of CuNCN with blue Cu, 
green N and gray C atoms. Exchange couplings are denoted 
by Ji, Ji and J', (b) Triangular arrangement of Cu atoms in 
the ab plane. 
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our local-probe study - including muon spin relaxation 
(/iSR), nuclear magnetic resonance (NMR) and electron 
spin resonance (ESR) - in the temperature range be- 
tween 63 mK and 300 K. Around 80 K an unconventional 
magnetically disordered phase develops implying extreme 
sensitivity of the magnetic state towards external pertur- 
bation. Our results bring forward both similarities and 
differences between the copper carbodiimide and oxides. 

We first turn to /iSR, which is an extremely sensitive 
probe for internal magnetic fields. It can easily distin- 
guish between fluctuating and static magnetism on one 
hand, as well as between LRO and static disorder on 
the other [lf|. In Fig. [2] we show /iSR results in zero 
magnetic field (ZF) and in weak transverse field (wTF) 
of 2 mT, measured on high-purity powder samples [17[ 
at the MUSR facility (ISIS, Rutherford Appleton Labo- 
ratory, UK) and at the GPS facility (Paul Scherrer In- 
stitute, Switzerland). Above 80 K, the ZF muon-spin 
polarization relaxes according to the relaxation function 
G(t) = Gvkt • exp[— At] in the Gaussian (j3 = 2) limit, 
with the Voigtian Kubo-Toyabe (VKT) function [H 
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VKT 



(*) 



+ -[l-(A^]exp[-(A^//3]. (1) 



The dominant relaxation is due to static nuclear- 
magnetic- fields distribution with the width A/7 M = 
0.16 mT ( 7 ,, = 2tt ■ 135.5 MHz/T), while A = 0.007 /is" 1 
represents a weak dynamical relaxation due to fast fluc- 
tuations of electronic magnetic fields. 

Below 80 K the shape of the ZF relaxation function be- 
comes irregular. The amplitude of the high-temperature 
component decreases at the expense of another, quickly 
relaxing component [Fig. (Ha)]. The latter prevails below 
20 K where the data again fit well to the VKT relaxation 
function [Eq. (JTJ], this time with f3 = 1.44 and signifi- 
cantly increased static internal fields, A/7 M = 17.7 mT 
[inset to Fig. Eta)]. The largely enhanced A can only 
be attributed to electronic fields. The observed dip and 
the ZF polarization leveling at 1/3 for longer times [l9[ 
affirm that local fields are disordered and fully frozen on 
the muon time scale [l6| . We stress that absence of co- 
herent oscillation of polarization discards the possibility 
of a LRO state down to at least 63 mK [inset to Fig. [21(a)] . 

The wTF measurements further disclose the magnetic 
state of CuNCN between 80 K and 20 K, where multi- 
component ZF relaxation is observed. The amplitude of 
the oscillating polarization, proportional to the fraction 
of muons still experiencing dynamical electronic fields, 
gradually decreases in this range. This witnesses pro- 
gressive development of static internal magnetic fields of 
electronic origin ["frozen" fraction in Fig. HJb)], whose 
magnitudes are large compared to the applied field. This 
phase is inhomogeneous as part of the muons detects 
static and the rest dynamical local fields. The wTF mea- 
surements confirm a fully static phase below 20 K. 

The single-dip ZF muon relaxation observed in 
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FIG. 2: (color online), (a) Relaxation of muon spin polariza- 
tion in zero magnetic field. Solid lines represent fits to Eq. JTJ, 
see text, (b) T-dependence of the "frozen" fraction, deduced 
from weak-transverse-field measurements. 



CuNCN below 20 K is commonly encountered in spin 
glasses (SG) 0. SG states have been detected by /xSR 
in several frustrated AFMs [201422T ] and appear due to 
some sort of randomness; structural disorder, random 
fields, bonds or strains, topological defects, etc. One 
could therefore hastily assign the CuNCN ground state 
as a SG state. Such state is generally characterized with 
a pronounced zero-field-coolcd/ficld-coolcd (ZFC/FC) ir- 
reversibility in Xb below the freezing temperature. We 
have, therefore, searched for the ZFC/FC irreversibility 
fl7| and for signs of a SG transition in the ac susceptibil- 
ity in various applied field between 0.8 mT and 5 T, but 
could not find any. Another peculiarity of CuNCN, un- 
characteristic for SGs, is the anomalous stability of the 
inhomogeneous region of coexisting static (low-T) and 
dynamic (high-T) phase [Fig. E£b)]. 

Although low-T yuSR results in CuNCN reveal static 
frozen fields, it is the absence of thermal remanence 
and well defined freezing temperature that rules out the 
canonical SG scenario. A similarly inhomogeneous phase 
has been reported in the frustrated kagome compound 
vesignieite and reasoned to be due to a vicinity of a quan- 
tum critical point 
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On the other hand, in few spin- 
gap systems partially frozen disordered states were unex- 
pectedly detected by /xSR [HI, [25j . They were attributed 
to a muon locally perturbing exchange pathways thus 
prohibiting singlet formation and liberating spin degrees 
of freedom in its vicinity 26|, [27J. For slow dynamics, a 
local-probe muon would not distinguish such a state from 
collective spin freezing. This is plausible in CuNCN, as 
muons are likely to stop near bridging N atoms of the 
NCN 2 ~ groups, similarly as protons are positioned in re- 
lated M(NCNH) 2 (M = Fe-Ni) [Hj]. 

To avoid any disturbance of the CuNCN structure by 
external probes we next decided to employ 14 N (7 = 1) 
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FIG. 3: (color online), (a) Comparison of Gaussian broad- 
ening a of 14 N NMR spectra and bulk susceptibility \b (in 
5 T). Inset: Simulation (red) of selected experimental (blue) 
spectra, (b) Spin-lattice relaxation rate 1/Ti divided by T. 
Inset: T-dependence of (i) the stretch exponent and (ii) 1/Ti. 
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FIG. 4: (color online), (a) Comparison of bulk susceptibility 
Xb and ESR spin-only susceptibility xesr, both measured in 
0.33 T. Inset: (i) Fits (red) of selected ESR spectra (green) 
and (ii) T-dependence of the 2D RVB order parameter, (b) 
T-dependence of the ESR line- width and the g-factor (inset). 



NMR measurements because intrinsic 1 N nuclei also di- 
rectly detect local magnetic fields. To obtain a reasonable 
signal-to-noise ratio the frequency-swept spectra were 
measured in high magnetic field of 9.4 T [29|. At 250 K 
the NMR spectrum (Fig. [3]) exhibits a characteristic 
quadrupolar powder line shape, which allows fixing the 
quadrupolar interaction parameters, vq = 0.79(2) MHz 
and r)Q = 0.06(2). In order to simulate the low-T spec- 
tra we assume a convolution with a T-dcpendent Gaus- 
sian local field distribution with the width o~(T) origi- 
nating from the hyperfine coupling between the nuclear 
and Cu 2+ electronic moments, a starts increasing dra- 
matically below ca. 100 K [Fig. EJa)] but fails to scale 
with the moderate increase of \b at low T [Fig. EJa)], 
thus ruling out the distribution of hyperfine couplings as 
the source of the spectral broadening. Moreover, spin- 
spin relaxation measurements (not shown) discard an in- 
crease in the relaxation rate as a possible origin of the 
low-T broadening, which is therefore inhomogeneous and 
due to a broad distribution of local susceptibilities {static 
electronic moments) at low T and high fields. 

The T dependence of the 14 N spin-lattice relaxation 
rate 1/Ti, measured with an inversion recovery tech- 
nique, reveals the development of electronic spin dynam- 
ics since the fast spin-lattice relaxation times of only a 
few ms [Fig. EJb)] directly measure the local-field fluctua- 
tions weight at the Larmor frequency. The 14 N magneti- 
zation relaxation curves were analyzed with the stretched 
exponential model M z (t) - M oc cxp[-(t/Ti) a ]. Al- 
though 1/Ti shows no obvious anomaly at 80 K, the 
stretch exponent a clearly starts to decrease below 80 K 
and reaches a value of 0.5 at lowest T [inset (i) to 
Fig. Efb)]. This evidences that the distribution of re- 
laxation rates widens below 80 K — another hallmark of 
a magnetically highly disordered state. Around 20 K a 



1/Ti maximum is observed [inset (ii) to Fig. (2Kb)], indi- 
cating freezing of spin dynamics below this temperature. 
The freezing is moderate, far from opening of a spin gap 
which would require a suppression of 1 /T\T, not observed 
experimentally [Fig.[3Jb)]. A ZF spin gap could however 
be closed by the large applied field of 9.4 T Q. 

The electronic properties of CuNCN were, therefore, 
also investigated at much lower fields (B ~ 0.33 T) with 
continuous-wave ESR measurements in X-band {vl = 
9.6 GHz), following both FC and ZFC temperature treat- 
ments. The ESR spectra [inset (i) to Fig. HI a)] comprise 
two overlapping components; a broad Lorentzian compo- 
nent with the RT line-width of 280(10) mT and a ca. 10- 
times narrower Gaussian component. The RT g- factor 
value g = 2.07 of the broad component is common for 
Cu 2+ in a square planar environment [3^]. The broad- 
component ESR intensity, corresponding to the spin-only 
susceptibility, was calibrated at RT with a couple of stan- 
dards, Xesr ~ 9(1) ■ 10~ 5 emu/mol. It is the same as the 
bulk susceptibility, which proves that the broad compo- 
nent represents the intrinsic ESR signal. The narrow 
component, on the other hand, is ca. 100-times less in- 
tense and exhibits a Curie-like T-dependence. The total 
fraction of impurity spins contributing to it is estimated 
as only 0.04%. The intrinsic ESR component exhibits a 
minimum in the line-width and a maximum in ^-factor 
around 130 K [Fig. SJb)]. Both indicate growth of spin 
correlations [30j , which appear particularly strong below 
~ 80 K. The featureless T-dependence of both parame- 
ters and no obvious ZFC/FC dissimilarities of ESR spec- 
tra speak against a conventional SG transition, in line 
with the lack of usual susceptibility fingerprints of such 
a state. Moreover, the intrinsic xesr rapidly diminishes 
below 80 K, thus diverging from \b [Fig. Ufa)]. 

The suppressed xesr implies opening of a spin gap 
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and thus allows for an alternative interpretation of the 
magnetic state in CuNCN. For the anisotropic triangu- 
lar lattice, Hayashi and Ogata proposed that in ZF two 
spin-singlet states are formed at different temperatures, 
depending on the J%/ J\ ratio [3l|. Below the upper crit- 
ical temperature 9 C = 3Ji/8 a quasi-lD gapless RVB 
phase sets in, while a 2D gapped RVB phase appears at 
the lower critical temperature Q' c = J2/8(l — 2J 2 /3Ji) 
[2]. We have further developed the analytical suscepti- 
bility expression [15( for the 2D RVB phase and arrive 
at 



Mb 



(l-ca)f(3V2J 2V ) 
3V2(Ji£-J2V) 



ca g( 3A /2J 2 7 ? ) 1 , (2) 



where a = J 2 ?7/ 'Ji£, with £, r\ as the ID, 2D RVB order 
parameters, respectively; / is a Fermi distribution func- 
tion and c = 1 + ( 1 — ^j- + \B\na) [32J . Assuming con- 
stant £ = I/v^tt [Hf in the 2D RVB phase we estimate 
the exchange constants Ji = 2300 K and J 2 = 540 K and 
establish the T-dependence of the 2D RVB order param- 
eter [inset (ii) to Fig.gfa)] from xesr(T) for 0' c = 80 K. 

We note that both RVB phases are characterized by 
complete spin-pairing, in agreement with no magnetic 
neutron scattering in polarized experiments The 
growing spin correlations observed in ESR are also in line 
with the 2D RVB phase below 6' c . This phase, however, 
appears to be very fragile, as muons can locally distort 
singlets and lead to static local states within the spin gap 
A s below 9' c . Moreover, as A s = 3\/2«/2?7 is comparable 
to the applied magnetic field in NMR, it is reasonable 
to assume that the lack of the spin gap and the grow- 
ing static fields observed in NMR may be field induced. 
This calls to mind the recently discovered emergence of 
inhomogcncous moments in extremely weak fields on an- 
other triangular lattice [33[ and points to the importance 
of frustration which amplifies the role of small pertur- 
bations in selecting between competing spin-liquid and 
more conventional frozen phases. 

In conclusion, CuNCN, the first nitrogen-based ana- 
logue of CuO, with a surprisingly low magnetic suscep- 
tibility, lacks magnetic ordering down to at least 63 mK. 
A magnetically highly disordered state has been wit- 
nessed below 80 K by complementary local probes. Al- 
ternatively to a SG-like frozen state featuring an un- 
precedented broadness of inhomogcncous static and dy- 
namic phase, no well-defined freezing temperature and 
no conventional susceptibility-fingerprint SG features, a 
dynamical spin liquid is proposed. The fragility of the 
electronic state to perturbations seems to be a peculiarity 
of carbodiimides as opposed to more conventional oxides. 
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SYNTHESIS AND CRYSTAL STRUCTURE 

All preparative steps were based on high purity 
cducts (99.9999% Cu, 99.999% HC1, highly crystalline 
cyanamide, milli-Q water) or on educts of analytical 
grade (concentrated NH3 solution, Na2SOs). After hav- 
ing synthesized the precursor Cu4(NCN) 2 NH 3 ,[l] it was 
oxidized over night in aqueous solution to yield X-ray 
pure CuNCN of excellent crystallinity [2], with notice- 
able Bragg peaks even at highest angles (Fig. [S]) . Neu- 
tron diffraction data [3] did not show any indication of 
defect formation on any of the atomic sites. 



a Quantum Design SQUID magnetometer. Zero-field- 
cooled (ZFC) measurements were performed after cooling 
the sample to the base temperature in zero field and field- 
cooled (FC) measurements after cooling in applied mag- 
netic field. In Fig. [5] we show measurements in various 
applied magnetic fields B in the range between 0.8 mT 
and 5 T. There is no obvious ZFC/FC splitting emerg- 
ing at low temperatures that would be characteristic of a 
canonical spin glass transition in any of the applied fields. 
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FIG. 5: (color online). Observed (red), calculated (black), 
and difference (blue) intensities of the X-ray Rietveld refine- 
ment of single-phase CuNCN using strictly monochromatized 
Cu-Kai radiation (A = 1.54056 A) with a flat-sample geom- 
etry and an imaging-plate detector. The vertical positions of 
the Bragg reflections are given in green. The residual values 
are 0.038 (R p ), 0.052 (i? wp ) and 0.035 (i?Bra gg ). 




100 200 
T(K) 



300 



100 200 
T(K) 



300 



FIG. 6: (color online). T-dependence of ZFC and FC bulk 
susceptibility in various applied magnetic fields. 



BULK MAGNETIC SUSCEPTIBILITY 



Magnetic susceptibility (\ = M/B) was measured on 
a 90 mg powder sample between 2 K and 300 K with 
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